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Non-diagonal open spin-1/2 XXZ quantum chains by separation of variables: 
Complete spectrum and matrix elements of some quasi-local operators 



G. Niccolfl 



Abstract The integrable quantum models, associated to the transfer matrices of the 6-vertex reflection algebra for spin 
1/2 representations, are studied in this paper. In the framework of Sklyanin's quantum separation of variables (SOV), 
we provide the complete characterization of the eigenvalues and eigenstates of the transfer matrix and the proof of the 
simplicity of the transfer matrix spectrum. Moreover, we use these integrable quantum models as further key examples 
for which to develop a method in the SOV framework to compute matrix elements of local operators. This method has 
been introduced first in (Q~) and then used also in [2|, it is based on the resolution of the quantum inverse problem (i.e. 
the reconstruction of all local operators in terms of the quantum separate variables) plus the computation of the action of 
separate covectors on separate vectors. In particular, for these integrable quantum models, which in the homogeneous 
limit reproduce the open spin-1/2 XXZ quantum chains with non-diagonal boundary conditions, we have obtained the 
SOV-reconstructions for a class of quasi-local operators and determinant formulae for the covector-vector actions. As 
consequence of these findings we provide one determinant formulae for the matrix elements of this class of reconstructed 
quasi-local operators on transfer matrix eigenstates. 
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1 Introduction 



In this paper we analyze the lattice quantum integrable models characterized in the quantum inverse scattering 
method (QISM) |[3l- lfl6l by (boundary) monodromy matrices which satisfy the reflection algebra |[T8l - ll26l w.r.t. 
the 6-vertex R-matrix solution of the Yang-Baxter equation. The prototypical elements in this class of quantum 
integrable models are the open XXZ spin- 1/2 quantum chains, reproduced under the homogeneous limit of the 
representations of this reflection algebra on the spin- 1/2 quantum chains. In the special representations corresponding 
to the open XXZ chain with diagonal boundary matrices^! the system has been largely analyzed in the framework 
of the algebraic Bethe ansatz (ABA) I3]-|4] with results going from the spectrum^ |[T8l - ll2"6l up to the correlation 
function^] |[29l - ll3Ti . where the Lyon group method^] li34l - B7l has been generalized to the reflection algebra case in 
the ABA framework. 

The situation is more complicated in the case of general non-diagonal boundary matrices; technical difficulties 
arise and the ABA method can be applied only limitedly to non-diagonal boundary matrices which satisfy special 
constrains which allow the definition of reference states by gauge transformation^j [51, 52]. Even in this class of 
constrained non-diagonal boundaries the use of algebraic Bethe ansatz is complicated from the fact that two sets of 
Bethe ansatz equations are generally required in order to have some numerical evidence of the completeness of the 
spectrum description ll53l . The other problem in the ABA framework is for the computation of scalar products which 
is missing as soon as non-diagonal boundary matrices are considered. Some analysis for solving this last problem 
has been addressed in the papers ||54] 1551 where the partition function for the dynamical diagonal open case was 
considered, see also 11561 . According to the standard techniques in the ABA framework, this is the first step towards 
the determination of scalar product formula in a determinant form. The situation for the general unconstrained non- 
diagonal boundary matrices is even more fragmented in the framework of Bethe ansatz analyst; only the eigenvalue 
analysis is implemented by the fusion procedure for special values (roots of unit) of the anisotropy parameter Il62l 1631 
while for the other cases (non-roots of unit) the construction of the Q-operator is mainly at a conjecture level [64 1. 
Once again the completeness of the spectrum description is verified only by some numerical analysis and the absence 
of eigenstates construction is the first fundamental missing step toward the matrix elements of local operators. 

The circumstances that for general non-diagonal boundary matrices the ABA method does not work while for the 
constrained ones, for which it works, it is instead missing a scalar product formula have so far prevented to further 
generalize the method described in [29]-[31 ] to the reflection algebra representations for non-diagonal boundary ma- 
trices. However, we are in the position to develop a different approach based on the Sklyanin's quantum separation 
of variables (SOV) method^] 11691 - 11711 . which can be used to get the exact characterization of their spectrum (eigen- 
values & eigenstates) and the computation of their form factors. Indeed, Sklyanin's SOV is a more efficient method 



2 Numerical 2x2 matrix solutions of the reflection equations, located at the ending points of the quantum chain. 
3 See 1 17] for the analysis by nested ABA of higher rank open spin chains with diagonal boundary conditions. 

4 Let us recall that in the half-infinite volume case (with one boundary) similar multiple integrals formulae have been previously derived 
in [ 27 1-| 28 1 by using the q- vertex operator method. 

5 Always in the ABA framework, see also [32 |-|33 ] for the extension of this method to the higher spin quantum chains. 

6 See the series of papers [48|-[50| for the analysis of XXZ spin 1/2 open chains with special constraints on the non-diagonal boundaries 
by a method combing coordinate Bethe ansatz and matrix ansatz methods. 

7 It is worth mentioning that by using a different approach, the representation theory of the so-called q-Onsager algebra, the spectrum of 
the spin- 1/2 open XXZ quantum chains with the general non-diagonal boundary conditions has been characterized in 11571 - 161 1 in terms of 
the roots of certain characteristic polynomials. 

8 Let us comment that in the special case of the spin- 1/2 representations of the rational 6-vertex reflection algebra the construction of the 
functional version of separation of variables of Sklyanin has been implemented in |65 66 1; that is a representation of the rational reflection 
algebra on a space of symmetric functions. However, the explicit construction of the SOV representation and of the transfer matrix eigenstates 
in the original Hilbert space of the quantum chain are not provided. See also [67 1 for the analysis by the functional SOV of the related but 
more general spin-boson model introduced and analyzed by ABA in 1681 . 
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to analyze the spectral problem w.r.t. other methods^ like the algebraic Bethe ansatz. It works for a large class of 
integrable quantum models; it leads to both the eigenvalues and the eigenstates of the transfer matrix providing a 
complete characterization of the spectrum under simple requirement^^ Moreover, in all the integrable quantum 
models analyzed by SOV in the series of papers [ 80 1-[ 85 T ] it was possible to show that the transfer matrix forms a 
complete^ set of commuting conserved charges of the quantum model. 

The approach used in the present article can be considered as the generalization to the SOV framework of the 
Lyon group method. It has been already implemented in [2| for the XXZ spin 1/2 quantum chain 11721 . [86]-[93] 
with antiperiodic boundary conditions. The results^ obtained in (21 go from the complete characterization of the 
spectrum up to the calculation of the form factors of the local spin a% in a determinant forrrj^l This approach has 
been originally introduced in |£TJ for the lattice quantum sine-Gordon model |@J [16] and then generalized in ||96l to 
the r 2 -modeO gTJ and the chiral Potts model ll98l - lfT09l . 

Finally, let us comment that the analysis of these systems with non-diagonal boundary conditions is in particular 
interesting for the relevant physical applications to systems in non-equilibrium like the asymmetric simple exclusion 
processes (ASEP) as they allow to describe systems for which the particle number is not conserved^!. Then, the lack 
of knowledge on the spectrum in this general framework, the complications emerging in the use of algebraic Bethe 
ansatz plus the physical relevance of these open chains with non-diagonal boundary conditions makes clear how 
big can be the impact brought to this research area from the solution of these systems by our approach in quantum 
separation of variables. 



2 Reflection algebra and open spin- 1/2 XXZ quantum chain 

In this section we describe a class of quantum integrable models characterized in the framework of the quantum 
inverse scattering method by monodromy matrices U(\) which are solutions of the following reflection equation: 

R 12 (A - fji) Ui (A) R 12 ( A + fj.) U 2 (/i) = U 2 (/i) Ri 2 (A + fi) U x (A) R 12 (\-(m), (2.1) 

w.r.t. the 6-vertex trigonometric solution of the Yang-Baxter equation, the i?-matrix: 

/ sinh(A + r/) \ 

^ _ sinhA sinh?? 

sinhr/ sinhA 

\ sinh(A + ??) / 

where R x . ~ C 2 is a 2-dimensional linear space. 

'Other important examples are the coordinate Bethe ansatz 0721 . 1731 and 1741 . the Baxter Q-operator method 0731 and the analytic Bethe 
ansatz f75l-f76l. 

l0 Note that on the contrary in the ABA framework a proof of completeness has been given only for some models; for example see 0771 
for the XXX Heisenberg model, [78 1 for the infinite quantum XXZ spin chain with domain wall boundary conditions, |79 | for the nonlinear 
quantum Schroedinger model and the references contained in these papers. 

"Note that this is the natural quantum analogue of the definition of classical complete integrability. 

12 Note that previous results on this model were the Baxter Q-operator [94] and the Sklyanin's functional separation of variables for the 
XXX chain |70) extended in |95 | to the XXZ case. 

13 The same type of results are derived in [84 | for rational 6-vertex Yang-Baxter algebra representations on antiperiodic spin-s quantum 
chains. 

l4 Note that in II 101 - 1 1 131 a previous analysis by SOV method of the T2 -model has been implemented. 

15 For relevant references on this subject from the point of view of the connection to quantum integrable models see [114|-[lT9j and 
reference there in. 



G End(Ri <8>R 2 ), (2.2) 
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2.1 H.w. representations of 6- vertex reflection algebra on spin-1/2 chains 

Let K(X; £, 5, r) be the following (general non-diagonal boundary) matrix: 

1 ( sinh(A + C) Ke T sinh2A \ 

K{XXAt) = ^C\ *r- S inh 2A sinh(C -A) J' (23) 

where £, k and r are arbitrary complex parameters, it is the most general scalar solution of the 6-vertex trigonometric 
reflection equation: 

R 12 (X - /i)ET 1 (A) J Ri 2 (A + n)K 2 (ji) = K 2 (fi)R 12 (X + f i)K 1 (X)R 12 {X - //). (2.4) 
Then following Sklyanin BUI , it is possible to construct in the 2 N -dimensional representation spaced: 

= ®n=l*n, (2-5) 



two classes of solutions to the same reflection equation ([2.4]) . In order to do so, let us define: 

K_(A)=ir(A-7//2;C-,<5_,T_), = K(X + r?/2; C+, <5+, r+), (2.6) 

where £±, <5±> t± are arbitrary complex parameters and the (bulk) monodromy matrix^: 

Af (A) = i? N (A - £n - rj/2) . . . R 02 (X - & - r,/2) R 01 (A - 6 - tj/2), M(X) = (-1) N a£ M*° (-A) a y , (2.7) 
Mo(A) G End(Ro <g> TZn), solution of the 6-vertex Yang-Baxter equation: 

i? 12 (A - /i)M 1 (A)M 2 (/i) = M 2 ( f j,)M 1 (X)R 12 (X - n). (2.8) 
Now we can define the following (boundary) monodromy matrices U±(X) € End(Ro ® 7£n) as it follows: 

WL(A) = M (A)iv-_(A)M (A)= ( £ ^ ^(A) J ' M 



««?(A) = <(A)J^(A)J(#(A)=( ^ 



(2.10) 



then U- (A) and V+(A) = UP (—A) define two classes of solutions of the reflection equation (12. 4p . 

As standard in the quantum inverse method and as it was proven in [20], from these monodromy matrices it is 
possible to define a commuting family of transfer matrices T(A) € End(7?fj) as it follows: 

T(X) = tr {K + (X) M(X) K-(X)M(X)} = ti {K + (X)U-(X)} = tr {if_(A)W + (A)}. (2.11) 

The problems that we address in this paper are the complete characterization of the spectrum (eigenvalue & eigen- 
states) of this transfer matrix and the computation of some matrix elements of quasi-local operators for two quite 
general classes of non-diagonal boundary matrices K±(X). It is then worth recalling that the open spin-1/2 XXZ 



"The representation space of a spin 1/2 quantum chain of N local sites each one associated to a 2-dimensional local space R n . 
l7 Note that here we have chosen a shifted definition of the inhomogeneity w.r.t. the one used in the articles 11291 1311 : in this case the 
homogeneous limit corresponds to f m = for m = 1, . . . , N. 
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quantum chain, with the most general non-diagonal integrable boundary conditions, is characterized by the following 
Hamiltonian: 

N-l 



H ND = + <? y i<T v i+l + coshf7ofaf +1 ) 



i=l 



Slllll 77 

[a'i cosh + 2k_ (erf cosh t_ + ia\ sinh r_ 



sinh £_ 



Slllll 77 

H r^[( a N coshC + + 2k_i_((Tm coshr + + icrK sinhr + ). (2.12) 

smh 

Hamiltonian which is reproduced in the homogeneous limit by the following derivative of the transfer matrix (|2.1ip : 

2(sinhr ? ) 1 - 2N d 
tr{if+ ( 77 /2)}tr{iir^(r ? /2)} dA ' v "'| A =,/a 



g NJ> = 3717+ 71 /.m Z7 ^- 71 /om — r ( A )| +constant. (2.13) 



2.2 First fundamental properties 

Here, some important properties about the generators of the reflection algebra A±(X), B±(X), C±(X) and T>±(\) are 
given as they will play a fundamental role in the solution of the transfer matrix T(A) spectral problem. 

Proposition 2.1 (]U- -reflection algebra). In the reflection algebra generated by the elements of U-{\) the quantum 
determinant: 

det£Y_(A) = sinh(2A-2r ? )[^4_(A + ??/2)^ + (-A + ? ? /2) -S_(A + r ? /2)C_(-A + ?7/2)] (2.14) 
i 

= sinh(2A-2r ? )[P_(A + 7 ? /2)P_(-A + 7?/2) - C_(A + r?/2)B_(-A + 7y/2)], (2.15) 

is central: 

[detW_(A),W_(//)] = 0. (2.16) 
i 

Moreover, it admits the following explicit expression: 

detW_(A) = sinh(2A - 2t?)A_(A + r//2)A_(-A + rj/2), (2.17) 
<? 

where: 

N 

A_(A) = <7_(A)a(A)d(-A), d(X) = a(X - V ) , a(A) = JJ sinh(A - £ n + rj/2) , (2.18) 

n=l 

and: 

_ sinh(A + a± ± r?/2) cosh(A + fa db rj/2) (2 J9) 

sinh a± cosh (3± 

where a± and j3± are defined in terms of the boundary parameters by: 

sinh a± cosh f3± = — — , cosh a± sinh f3± = — . (2.20) 

2k± 2k± 

Moreover, the generator families „4_(A) and "D_(A) are related by the following parity relation: 

^ sinh(2A — ri) . , _ sinh 7? . ... 

while for the other two families the following parity relations hold: 

sinh(2A + 77) . \\_ sinh(2A + r/) 

^- (_A) " "sinh(2A-»;) B - (A) ' C " (_A) " "sinh(2A - rff'^ ^ 
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Proof. This proposition is just a rephrasing and a simple extension to the case of general non-diagonal if- (A) 
boundary matrix of the results stated in Propositions 5, 6 and 7 of Sklyanin's article [20]. The quantum determinant 
as defined in formulae (38)[2Qj-(42)[2Q|: 

detZY_(A) = U-(X-r]/2)U-(\ + r]/2) = W_(A + t?/2)ZL(A - 77/2) (2.23) 
i 

= A-(\ + r]/2)V-(\-ri/2) -B-(\ + ri/2)V-(\-r]/2) (2.24) 
is central independently from K- (A) being diagonal or non-diagonal. Here, we have used the definition (43) [20]: 

U{x) = tr 12 R 12 (- v ) (U-) 2 (X)R 21 (2X) = ( P_(A) -B_(A) \ (2 ^ 



sinh 77 ^ -C_(A) ,4_(A) 

P_ (A) sinh 2A - A- (A) sinh ry - sinh(2A + r?)fi_ (A) 

- sinh(2A + rj)B- (A) „4_ (A) sinh 2A - P_ (A) sinh 77 



So, we can also write it as it follows: 



C(A) Z>(A) 



and 



(2.26) 



detU-{\) = detK_(A)detM (A)detM (-A), (2.27) 
9 <?<?<? 

where det g M(A) is the (bulk) quantum determinant of the Yang-Baxter algebra: 

det M (A) = A(X + n/2)D(X - ri/2) - B(\ + rj/2)C(X - n/2) 

= a(X + r]/2)d(X-r)/2). (2.28) 

Here, we have denoted: 

Mo(A)J^ *£>Y (2.29) 



det tf_(A) = (A + 77 /2)( J fiC_) 2)2 (A - n/2) - (tf_) 12 (A + r,/2)(K„) 2A (X - t?/2) 

= - Smh(2A ~ 27? ) (sinh(A + C-) sinh(A - C-) + n 2 _ sinh 2 2A), (2.30) 
sinh £_ 

where: 

Kl(X) = sinh(2A - r/^ivl (-A)o#. (2.31) 
Then the explicit expression (12.17t follows observing that it holds: 

detif_(A) = sinh(2A-2r/) 5 _(A + r//2) 5 _(-A + 7 ? /2), (2.32) 
<? 

when we use the parameters a_ and /?_. 

Finally, it is simple to remark that formula (46) [20 is equivalent to the symmetry properties (I2.2U and (12.22 \ , 
which in turn imply the expressions for the quantum determinant (12.141 ) and (12- 15b . when used to rewrite formula 
(12241 . □ 

It is worth remarking that similar statements hold for the reflection algebra generated by U+(X). In fact, they are 
simply consequences of the previous proposition when it is taken into account that U^(— A) satisfies the same 
reflection equation ofU-(X). 
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Proposition 2.2 (£Y + -reflection algebra). In the reflection algebra generated by the elements of £/+(A) the quantum 
determinant: 

&etU+(X) = sinh(2A + 2rj)[A+(X - r]/2)A + (-\ - rj/2) - B+{X - rj/2)C + (-X - r//2)l (2.33) 
i 

= sinh(2A + 2 V )[V + (-X - r]/2)V + (X - V /2) - C+(-X - r 1 /2)B + {\ - rj/2)], (2.34) 

is central: 

[det« + (A),Z4(V)l =0. (2.35) 

i 

Moreover, it admits the following explicit expression: 

det^+(A) = sinh(2A + 2r ? )D+(A - r//2)D + (-A - rj/2), (2.36) 

w/jere the function D + (A) is defined by: 

D+(X) = g+(X)a(-X)d(X), (2.37) 

where g+(X) is defined in ()2. 19j) . Moreover, the generator families .4+ (A) and T> + (X) are related by the following 
parity relation: 

^^^^^-^-^ 

while for the other two families the following parity relations hold: 

8 +<- A > = - Mlml'V w • c +<- A > = - Si "!:^I"! c + (»- 0.39) 

smn(2A + rj) smn(2A + r\) 

Proof. Here, it is only need to notice that: 

det K + (X) = (K + ) hl (A - r ? /2)(i^) 2 , 2 (A + V /2) - (K + ) 1>2 (A - 7? /2)(Z+) 2)1 (A + r?/2) 

sinh(2A + 2r/) 
sinh 2 

where: 



-(sinh(A + C+) sinh(A - ( + ) + ^ sinh 2 2A), (2.40) 



K+(X) = sinh(2A + V )a y K^(-X)a y , (2.41) 

can be written in the form: 

detK + (X) = sinh(2A + 2r])g + (X - r}/2)g + (-X - rj/2)), (2.42) 

where the parameters a + and /3 + , entering in the function g + (A), are defined in terms of the parameters of the K + (A) 
boundary matrix in (j2.20p . □ 

Let us introduce the following notations: 

K m _ 1 ( sinh(A + C± ± r?/2) K±e T± sinh(2A ±rj) \ _ f a± (A) 6± (A) \ 
±l J sinhC± I K±e- r ± sinh(2A ± r/) sinh(C± T r?/2 - A) i I c± (A) d± (A) i ' 1 ' J 



and by using these notations let us rewrite the transfer matrix (|2.1ip in the following two equivalent forms: 



T(X) = 77 (±) (A) + b T (X)C ± (X)+c T (X)B±(X), (2.44) 



where: 

r W( A ) = aT ( A ) A± ( X ) + d T (A) V±(\), (2.45) 
are the transfer matrix0of the system with diagonal matrix K T (X), respectively, then: 
Corollary 2.1. 7^ (A) admits the following explicitly even forms w.r.t. the spectral parameter X: 



T\ (±) (A) = a T (A)A t (A) + a T (-A)A t (-A) (2.46) 
= d T (A)P±(A)+d T (-A)P±(-A) ; (2.47) 



a±(A) ee sinh ( 2A ± V) sinh ( A + C± ± rj/2) 
sinh 2A sinh C,± 

^ S inh(2A±7/) S inh(C±-A±77/2) ^ (2 49) 
sinh 2A sinh 

Moreover, also the most general transfer matrix is even in the spectral parameter A: 

T(-A) = T(A). (2.50) 

Proof. By using the formulae (|2.2ip and (|2.38p to rewrite 7y (A) only in terms of A±(\) or only in terms of 
T>±(\) after some simple algebra we get our formulae (|2.46p and (|2.47p . respectively. Then the parity (|2.50p of the 
transfer matrix T(A) follows remarking that the parity properties: 

b T (-\)C ± (-\) = b T (\)C±(\), c T (-\)B±(-\) = c T (\)B ± (\), (2.51) 

are just a rewriting of the known properties (|2.22|) and (I2.39p . □ 

Proposition 2.3. The monodromy matrix U±{\) satisfy the following transformation properties under Hermitian 
conjugation: 

I) Under the condition rj € iM. (massless regime), it holds: 

W±(A)t = [W±(-A*)]*°, (2.52) 

for {» T±J iK±, iC±, £i, Cn} G K n + 3 . 

II) Under the condition r) £ K (massive regime), it holds: 

U±(\)^ = [U±(\*)] t0 , (2-53) 

/or {r±,K±,C±,*6, GM N+3 . 
Under the same conditions on the parameters of the representation it holds: 

T(A) t = T(A*) ; (2.54) 

z'.e 7~(A) defines a one-parameter family of normal operators which are self-adjoint both for A real and imaginary. 



I8 Note that 7^ + '(A) corresponds to _K_(A) diagonal while K+(X) is left general as well as 77 '(A) corresponds to K+(X) diagonal 
while K- (A) is left general. 



9 



Proof. This proposition gives the generalization to the case of general non-diagonal boundary conditions of the 
transformation properties under Hermitian conjugation proven in [29] for the diagonal case. The proof is given by 
using the following transformation properties under Hermitian conjugation: 

I) For {irj, ir±,iK±, «C±,6,---,64 S M N+3 then it holds: 



flon(A-£i-»7/2)t 



cr y Ron(X* ~ 6 - r]/2)a y = - [i? „(-A* + £1 - r,/2)] to 
[K ± (-\*)] t0 ; 



II) For{77,r±,K±,C±,^i) ...,^n} S M N+3 then it holds: 



flon(A-6-r?/2) f 



-ffglWA* - 6 - »7/2)ag = [i? n(A* + 6 - r?/2)]*° 
[tf±(A*)]<°; 



(2.55) 
(2.56) 



(2.57) 
(2.58) 



which can be verified by direct calculations. From these it follows: 

N 



v_ 

if 



M(A)t = 

for the (bulk) monodromy matrix and so: 



7]' 



\ ) a*K 



o 



^±(A) t 



M( 



to 



to 



M( 



to 



to 



q-reverse-order 



(2.59) 



(2.60) 



where the notation q-reverse-order is referred to the reverse order in the generators of the Yang-Baxter algebras. 
More in details, the matrix elements of the 2x2 matrix li±{\y are computed by using the normal matrix products 
in the auxiliary space 0, as indicated inside the brackets at the r.h.s. of (12.601 ). then in each element of the matrix 



to 



to verify that the r.h.s. of (I2.60l l with this prescribed order coincides with: 



to 



W±(A)t the matrix elements of M( ( ) A*) are put to the right of those of M( ( ) A*) . It is then simple 



(2.61) 



"=U^)A J > 



to 



which proves both (12.521 ) and (12.53b - Then by using these last two formulae and the transformation properties (12.56) 
and (12.581 ). we get: 



T(A)t = tr 



to 



T((—)\*) = T(A* 



(2.62) 
□ 



3 SOV representations for T( A) -spectral problem 

The method to construct quantum separation of variable (SOV) representations for the spectral problem of the trans- 
fer matrices associated to the representations of the Yang-Baxter algebra has been defined by Sklyanin in 116911701 im . 
Here, we show that there are quite general representations of the reflection algebra with non-diagonal boundary ma- 
trices for which the quantum SOV representations can be constructed by adapting Sklyanin's method. More in details 
the following theorem holds: 
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Theorem 3.1. Let the inhomogeneities {£1, £n} € C satisfy the following conditions: 

ia^ib + rr] Va^G{l,..,N}Wr€{-l,0,l}, (3.1) 

then the commuting families of generators of the reflection algebra B±(\) and C±(X) are diagonalizable and with 
simple spectrum. 

Moreover, denoted with e € {+, — } the following statements hold: 

I) The representations for which the commuting family B t {\) is diagonal define the quantum SOV representations 
for the spectral problem of the transfer matrix T(A) associated to the boundary matrix K_ e (A) diagonal or lower 
triangular while K t (A) is a general non-lower triangular; i.e. the SOV representations for: 

T(A) = T^\\) + c_ e (A) B e (A), for 6_ £ (A) = and b e (A) ^ 0. (3.2) 

II) The representations for which the commuting family C e (\) is diagonal define the quantum SOV representations 
for the spectral problem of the transfer matrix T(A) associated to the boundary matrix K^ t (A) diagonal or upper 
triangular while K t (A) is a general non-upper triangular one; i.e. the SOV representations for: 

T(A) = 7^ {e) (A) + 6_ e (A) C e (A), for c_ e (A) = and c e (A) + 0. (3.3) 

The proof of the theorem will be given in the following sections by the explicit construction of £>±(A)-eigenbasis and 
the solution by quantum separation of variables of the spectral problem for the transfer matrix T(A) in the class of 
representations of the reflection algebra listed in point I) of the theorem. The Hermitian conjugation properties lead 
to the construction of the C± (A)-eigenbasis and then imply that the theorem holds also for the class of representations 
listed in point II) of the theorem. 



3.1 Left and right representations of the reflection algebras 

Let us give some more details on the space of the representation of the spin 1/2 quantum chain. As this is the same 
representation space used in [2 ] for the 6-vertex Yang Baxter algebra similar notation will be given here. 

Let us introduce the standard spin basis for the 2-dimensional linear space R n , the quantum space in the site n of the 
chain, whose elements are the a z n -eigenvectors \k, n), characterized by: 

a z n \k,n) = k\k,n), fee {-1,1}. (3.4) 

Similarly, the cr,^-eigenco vectors {k, n\, characterized by: 

<Tn\k,n) = k\k,n), fee {-1,1}, (3.5) 

define a basis in L n the dual space of R n . Then, 2 N -dimensional representations with N + 6 parameters (the inhomo- 
geneities and the boundary parameters) of the reflection algebra are defined in the left (covectors) and right (vectors) 
linear spaces: 

£n = ®n=iL„, TZ N = ^ =1 R n - (3.6) 

Moreover, TZn is naturally provided with the structure of Hilbert spece by introducing the scalar product character- 
ized by the following action on the spin basis: 

N 

(®" =1 \k n ,n),^ =1 \k' n ,n)) = ]j5 hnjK Vk n ,k' n e {-1,1}. (3.7) 

n=l 
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3.2 B -SOV representations of reflection algebra 



In this subsection we construct the left and right SOV-representations of reflection algebra generated by U-{\) by 
constructing the left and right £>_(A)-eigenbasis. 

Theorem 3.2. I) Left B-(X) SOV-representations If (|3,ip is satisfied and 6_ (A) 7^ 0, then the states: 

N 



h hi- 1 tn\T\ ( A -W 2 -M \ hn 



where: 



(0|=^ =1 (l,n|, N_ 



(i) Jih 



ril<fe<a<N( r /a Va 



and 



(0| nn=i^-(^/2-e.)/A-(r7/2-Cn) |0) 



1, 



1/2 



(£n + (An - g)^ 



== cosh 2 

h n G {0, 1}, n G {1, N}, define a B-{X)-eigenbasis of C^: 

<-,*|5_(A) = b_,»(A)(-,* 
where (—, h\ = (—, hi, h^\forh= (hi, /in) and 

T _ sinh(2A — rj) 



B_ A (A) = K_e T 



sinh 



a/,(A)a/,(-A), 



a/, 



_ 1 
(A) = J] sinh(A - £ n - (h n - -)r/). 



n=l 



f?« f/je generic state (—, h\, the action of the remaining reflection algebra generators follows by: 



mi^-s S^ n 



cosh 2A — cosh 2C 



C»6) 



^ sinh(2Cr° j - r?) sinh 2^" a; 6 =1 cosh 2C«1 M - cosh 2Q 

by^a modN 

, M , ^ , „ , x - r \cosh2A-cosh2cf b) , 
x (-,*|77*'« + detM(O) cosh(A - 77/2) J] 

9 6=1 cosh 77 — cosh 2Q 

/ xN / ,x , x - r \cosh2A-cosh2d' li,) , 

+ (-1) N cothC- det M(m/2) sinh(A - 77/2) TT lirr(- * 

« 6=1 cosh 77 + cosh 2Q bj 



h a (rt h °h 



where: 



<Pr, 



£n + (K - -)r] 



for h n G {0, 1} and Vn G {1, 2N}, 



1 - 26(a - N) with 6(x) = {Ofor x<0, I for x > 0}, 



and: 



(—,hi,...,h a ,...,hfi\T^ 



, hi, h a ± 1, h^i I 



(3.8) 



(3.9) 



(3.10) 



(3-11) 



(3.12) 



(3.13) 



(3-14) 

(3.15) 
(3.16) 

(3.17) 



Indeed, the representation of V '-(A) follows from the identity \2.2l \ while C_(A) is uniquely defined by the quantum 
determinant relation. 
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II) Right B-(\) SOV-representations If (|3. 1 j) is satisfied and 6_ (A) 7^ 0, the states: 

(i-M 



where: 



h n G {0, 1}, n G {1, N}, define a B^(\)-eigenbasis of TZ^: 

B_(A)|- A) = |-,A)B_, Jk (A). (3.20) 
On f/je generic state \— , h), the action of the remaining reflection algebra generators follows by: 

S sinh(2Cr a) -r ? )sinh2Cr a) b= i cosh 2^ - cosh 2^ 

, x Acosh2A-cosh2d M 
+ \-,h) det M(0) cosh(A - rj/2) ]J 7^7- 

q b=1 cosh 77 — cosh 2Q b 

+ (_i)N+i|_ /,) coth c _ det Af (»tt/2) sinh(A - v /2) 17 cosh 2A ~ cosh 2 Cf ^ (3 21) 

9 b=l cosh r/ + cosh 2Q bJ 

D-(Cf a) ) = (A;H)^A_(C^ - 2^a), Tt\-,h u h a , h N ) = [-, h x , h a ± 1, fc N ). (3.22) 

Indeed, the representation of A- (A) follows from the identity \2.21\ while C_(A) w uniquely defined by the quantum 
determinant relation. 

Proof of I). It is worth writing explicitly the (boundary-bulk) decomposition of the reflection algebra generator: 

B_(A) = -o_(A)A(A)S(-A) + 6_(A)A(A)A(-A) - c_(\)B(\)B(-\) + d_(A)B(A)A(-A), (3.23) 
in terms of the generators of the Yang-Baxter algebra. Then, the following well know properties: 

(0|A(A) = a(A)(0|, (0|D(A)=d(A)(0|, (0|fl(A) = 0, <0|C(A) + 0, (3.24) 

with: 

N N 

o(A) = [Jsinh(A-^ n + 7//2), d(A) = J] sinh(A - £ n - T//2) , (3.25) 

n=l n=l 

imply that (0| is a 0_(A)-eigenstate with non-zero eigenvalue: 

(0|£_(A) = b_ j0 (A)<0| (3-26) 

where: 



B (A) = K_e- sinh ( 2A 7 ?) Q ( A )a(-A). (3.27) 
smh C- 
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Now by using the reflection algebra commutation relations: 

A-MB-M = Smh(Al " " 2 \"\ S[n ln 2 + i\~ ^ B-i^A-^) 
' v ' smh(Ai - A 2 )smh(Ai + A 2 ) 

sinh(2Ai -77) sinh 77 

+ S inh(A 2 -A 1 )sinh2A 1 ^ (A2) " 4 - (Al) 

SUlh// B_(A 2 )P_(A 1 ) (3.28) 



sinh(Ai + A 2 ) sinh 2Ai 

we can follow step by step the proof given in d to prove the validity of (|3.1ip and (|3.20p . Under the condition 
(|3,ip . these relations also imply that each set of states (— , h| and |— , h) form a set of 2 N independent states, i.e. a 
/3_(A)-eigenbasis of £n and IZ^, respectively. 

The action of A-(cj ) hb ^) for b € {1, 2N} follows by the definition of the states (— , h|, the reflection algebra 
commutation relations (I3.28P and the quantum determinant relations. Moreover, by using the identities: 

U-(r}/2) = detAf(O) I , U-(rj/2 + iir/2) = icothC- detM(ivr/2) oL (3.29) 
and remarking that A-(\) has the following functional dependence w.r.t. A: 

2N+1 

^_(A) = e^ a - m+1 ^A^ a (3.30) 

a=0 

we get the following interpolation formula for the action on (— , h|: 



2N 2N+2 • i q Jh b h 

( -,hi^w=E n . L S> )/ -(^')(-.i.iT a - 



-fa 



a=l h 



\ sinh(cf a) -Cf 6) ) 



sinh (A — Q 



+ detM(0) TT —--^ (_, h | 



6=1 sinh(Q N ; +1 - Q 

6^2N+l 



--- . - sinh ( A /"v" u '' ) 

+ icothC-detM(i7r/2) TT ^ — b ' (-,hl, (3.31) 

il 1 smh(cS + 2-(f b) ) 



where: 



2N+l = T 7A CS+2 = + W2, detM(A) = a(A + ? ? /2) a ((A-77/2), (3.32) 



and we have denoted Cm+b^ = ^2N+b f° r ^ = 1)2 for any £>_-eigenstate. Then, it is a simple exercise to rewrite 
this in the form (pTL4]) . □ 

Proof of II). The proof is given along the same line delineated for the point I) of the theorem, we just need to make 
the following remarks. First of all being: 

A(A)jO) = d(A)jO), £>(A)|0) = a(A)|0), S(A)|0)=0, C(A)]0) ^ 0, (3.33) 

then |0) is a £>_(A)-eigenstate with non-zero eigenvalue: 

fi_(A)M = B_, 1 (A)M (3-34) 

where: 

B !(A) = K„e T " smh (2A-^) (3 35) 

smh 
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Now all we need are the following reflection algebra commutation relations: 

sinh(Ai — A 2 + rj) sinh(A 2 + \\ - rf) 



B.{\i)V^{\ 2 ) 



sinh(Ai — A 2 ) sinh(Ai + A 2 
sinfi7/sinh(A 2 + Ai — rf) 



V_{\ 2 )B-{\i 



sinh(Ai — A 2 ) sinh(A 2 + Ai) 
sinh?7 



Z?_(Ai)B_(A 2 ) 



. , , , ,-^_(Ai)S_(A 2 ). 
smh(Ai + A 2 ) 

By using them we get the following interpolation formula for the action on |— , h): 



2N 



2N+2 



a=l 



+ |-,h)detM(0) Y\ 



b Ji shMd ha) - C ( b hb) ) 

sinh(A - d 



6=1 smh (C 2 N+l ~ Cfe 
b^2N+l 

2N+1 



(/*)> 



— 1\ 



T-r sinh(A - d hh) ) 
,h)cothC-detM^/2) J] . u l " b J bW 
i b=1 smh(^ N ; +2 -Q ') 



which can be rewritten in the form (|3.2ip . 



(3.36) 



(3.37) 



□ 



3.3 /3+-SOV representations of reflection algebra 

In this subsection we construct the left and right SOV-representations of the reflection algebra generated by U + (X) 
by constructing the left and right ,6 + (A)-eigenbasis. 

Theorem 3.3. I) Left B + (X) SOV-representations If (|3.ip is satisfied and b + (A) ^ 0, then the states: 



1 N (r> ( r (1 h\ (1 ~ M 

(+,h u ..,h K \ = ^(°\U l ] 



+ ^iVD+(-CA- 



where 



N+ 



ni 



<b<a<N 



(riF-rfl 1 *) 



1/2 



.<oi (nn=i^ + (-d 1) )/D + (-d 1) )) io> 

h n € {0, 1}, n G {1, N}, define a B + {\)-eigenbasis of C^\: 

(+,A|B + (A) = b +j »(A)<+,*| 1 
where (+, h\ = (+, hi, h^\for h= (hi, h^) and 

sinh(2A + 77) 

B+,*(A) = K + e T + ^— -a h {\)a h {-\). 

sinn Q + 



(3.38) 



(3.39) 



(3.40) 



(3.41) 
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On the generic state (+, h\, the action of the remaining reflection algebra generators follows by: 

-{h b ) 



^ sinh(2Cr a) + rj) sinh 2^ £=i cosh 2 G - cosh 2Cr j 



N cosh 2A -cosh 2& hb) 



(+,h\T£ a + detM(O) cosh(A + t//2) TT tttx+^I 

9 6 A = i cosh ?? - cosh ^ 



+ (-1) N+1 cothC+ det M(ivr/2) sinh(A + rj/2) TT cosh2A cosh 2 f " {+ ^ (342) 
9 b=i cosh r/ + cosh 2C 6 ( bj 

...,/i a , ...,h N \lf = {+,hi, ...,h a ± 1, ...,h N \. (3.43) 

Indeed, the representation of A+{X) follows from the identity < \2. 38i while C+(A) « uniquely defined by the quantum 
determinant relation. 

II) Right £>+(A) SOV-representations If w satisfied and 6+ (A) 7^ 0, f/zera the states: 

1 N / a rr (0 h \ hn — . 

|0> * 

_ sinh(2{ ,) 
smh(2£ n + r;) 

/i n G {0, 1}, n G {1, N}, define a B + {X)-eigenbasis oflZ^: 

B+(X)\+,h) = \+,h)B +th (X). (3.46) 
On the generic state |+, h), the action of the remaining reflection algebra generators follows by: 

^ +W i + .,>^^i + .'-) si " hpA ;f il ''' (A+ ^ n ^r^v, ^'' 

S sinh(2Cr a) +r ? )sinh2Cr a) £J cosh 2 Q ~ cosh 2 Q 

fe^a modN 

, N tt cosh 2A - cosh 2(9 hh ^ > 
+ \+,h) det M O cosh(A + V /2) TT ^ 

* 6= 1 cosh 77- cosh 2Cp j 

N cosh 2A — cosh 2C^ hb ^ 

+ (-l) N+1 |+>^> cothC+ det M(m/2) sinh(A + r?/2) TT ^— , (3.47) 

q cosh r? + cosh 2^ ft) 

A+(Ci M ) = D + (C^ a) - 2 ¥ > a &), 7^|+, /» N ) = |+, /n, fc a ± 1, h N ). (3.48) 

Indeed, the representation of T> + {\) follows from the identity (12. 38i while C+(A) is uniquely defined by the quantum 
determinant relation. 

Proof of I). The proof is given along the same line delineated in the previous theorem, we just need to do the follow- 
ing remarks. First of all let us give the (boundarydoulk) decomposition of the reflection algebra generator: 

B+(X) = B{X)D(-X)a+(X) + D(X)D(-X)b+{X) - B{X)B(-X)c + (X) - D{X)B{-X)d+(X), (3.49) 
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in terms of the generators of the Yang-Baxter algebra. Then, the properties (I3.24t imply that (0| is a £? + (A)-eigenstate 
with non-zero eigenvalue: 

B +)1 (A) ee s j^±Hl d{X ) d (-X). (3.50) 

The proof of the point I) is based on the following reflection algebra commutation relations: 

V + (M)B + (h) = ^MA,-A 2+ ,,)dnh(A 2 + A 1+ ,,) 

smn(Ai — A2) smh(Ai + A2) 
sinh(2A 2 + r?)sinhr/ 



sinh(Ai — A2)sinh2A2 
sinh?7sinh(2A2 + rj) 



sinh(Ai + A2) sinh 2A2 

and on the identities: 



B + (X 1 )V + (-X 2 ), (3.51) 



U+(-rj/2) = detM(O) I , U + {-ri/2 + iir/2) = i coth C+ det M(m/2) u%. (3.52) 
By using them and the fact that T> + (A) has the following functional dependence w.r.t. A: 

2N+1 

V+(X) = e (2a ~ 2N+ V x V + , a (3.53) 

a=0 

we get the following interpolation formula for the action on (+, h|: 

2N 2N+2 ■ u(\ _ Ahb)\ 

2N+2 . Jh b )s 



b Ji sinh(C«i» +1 - Cf" 1 } 
fe^2N+l 

2N+1 sinhCA - C (/lb)N ) 
! coth C+ det M(iw/2) T[ ) ^ b I (+,h|, (3.54) 



where we have denoted C2N+6(^2N+&) = C2N+& ^ or b = l,2 and 

cS+i = -^a cfita = -^/ 2 + i7r / 2 - < 3 - 55 ) 

Then, it is a simple exercise to rewrite this in the form (|3.42p . □ 

Proof of II). Similarly, by using (13.491 ) and (I3.33I ). it follows that |0) is a Z3+(A)-eigenstate with non-zero eigenvalue: 

B + (A) = K+ e r+ Smh ^ A , + T?) a(A)a(-A). (3.56) 
sinh <, + 

Then the proof of the point II) is based on the following reflection algebra commutation relations: 

1MA.VWA0 = si "MA,-A 2 + ,) S mh(A 2 + A 1 + ,) 

smh(Ai — A2J smh(Ai + A2) 
sinhT ? sinh(2A 2 + y) , , , 
smn(A2 — Ai)smh2A2 
+ B mh, S mh(2A 2 + ,,) A 
smh(Ai + A2) smh 2A2 
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By using them we get the following interpolation formula for the action on |+, h): 



2N ^i\i-r^ . , / , /-(' l bJ\ 

^l + (A)| + ,h)=^TM + ,h) n ,) b (H b )M & a) ) 

ti 6=1 smHd ha) - crn 

2N+2 . Jh,), 

+ |+,h)detM(0) TT m - 1 i 

b^2N+l 

2N+1 ginhfA _ Ah b ) ) 

+ i|+,h)cothC+detM(i7r/2) IT ^ ' , (3.58) 

£J sinh(c5 1 N ) +2 -C^ ) ) 

which can be rewritten in the form (|3.47p . □ 

4 SOV-decomposition of the identity 

The action of a generic left S-t-eigenstate on a generic right S±-eigenstate are here compute in this way allowing 
to write the decomposition of the identity in the corresponding basis. It is worth remarking that for the Hermitian 
conjugation properties of Proposition 12.3 1 these results correspond to the computations of scalar products between 
B± -eigenvectors and C± -eigenvectors. We show that up to an overall constant these are completely fixed by the left 
and right SOV-representations of the Yang-Baxter algebras when the gauge in the SOV-representations are chosen. 

4.1 Change of basis properties 

Let us present the main properties of the 2 N x 2 N matrices U^ L ' e ' and JJ^ R '^: 

2 N 2 N 

( e ,h| = (h|C/(^) = ^^(x- 1 «l and \e,h) = U^\h) = ^U i i R ^- 1 (i)), (4.1) 
i=l i=i 

which define the change of basis to the SOV-basis starting from the original spin basis: 

(h| = ®" =1 (2h n - l,n| and |h) = ®^ =1 \2h n - 1, n), (4.2) 
where xis the following natural isomorphism between the sets {0, 1} N and {1, 2 N }: 

N 

x:he{0,l} N 4/(h)El + ^ 2^h a G {1, 2 N }. (4.3) 

a=l 

Note that the matrices JJ^^ and U^ R ' e ^ are invertible matrices for the diagonalizability of B e {\): 

U {L ^B e {\) = A Bt {X)U {L ^\ B e {\)U {R ^ = U (R ^A Be (X). (4.4) 
Here Ag e (A) is the 2 N x 2 N diagonal matrix whose of elements, for the simplicity of the £> e -spectrum, read: 

(A Se (A)). ■ = ^B ejX _ 1(i) (A) Vi,j G {1,...,2 N }. (4.5) 

Moreover, it holds: 
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Proposition 4.1. The 2 N x 2 N matrix: 

M (e) = jjiL^jji^e) (4 _ 6) 

is diagonal and it is characterized by: 

l<b<a<N 'la '/{, 

Proof. Note that being the action of a left Z^-eigenstate on a right B e -eigenstate, corresponding to different B € - 
eigenvalues, zero this implies that the matrix is diagonal; then to compute its diagonal elements we compute 
the matrix elements 9^ = (— , hi, h a = 0, /in \A-(£, a +i]/2)\ — , hi, h a = 1, /in), where a € {1, N}. 
Using the left action of the operator A-(£ a + we § et: 

flH - a (n/2 n Sinhr? A cos^cP-cosb^C^ 
W Uj ^(26,-^)11 cosh 2C (0)_ cosh 2C (M 

x (-,hi,...,h a = l,...,h N \-,h 1 ,...,h a = l,...,/i N ) (4.8) 
while using the decomposition (12.211) and the fact that: 

V-i-ta -v/2)\-,hi,...,h a = l,...,h N }=0 (4.9) 

it holds: 

A-(Sa + T]/2)\-,h 1 ,-,ha = l,...M = ■ fe^TX N A -(W2~ga) 

smh(2f a +r/J 

x |-,/i!,...,/i a = 0,...,/t N ), (4.10) 

and then we get: 

°a~ )= ■ ul^ V ^ -(v/2-L)(-,hi,...,K = 0,...,h N \-,h l ,...,h a =0,...,h N ). (4.11) 
sinh(2£ a - rj) 



so that it holds: 



(-,hi,...,h a = l,...,h N \-,hi,...,h a = l,...,h N ) _ A cosh 2(i 0) - cosh 2(j y hb) 
(-,hi,...,h a = 0,...,h N \-,hi,...,h a = 0,...,h N ) ^ cosh 2^ 1} - cosh 2cf 6) ' 



from which one can prove: 



(-,hi,...,h N \-,hi,...,h N ) -j-r ija ] - 

<— , 1, 1, 1> 11 Jha)_(h b )- V- L *> 

1 K6<o<N Vb 



This prove the proposition for e = — , being 



<-, i, ii-, i, i> = n (1) 1 (4.i4) 

l<6<a<N % 

by our definition of the normalization N_. Similarly for e = +, we compute the matrix elements 9i + ^ = 
(+,hi,...,h a = 1, ...,h N \V + (€ a - rj/2)\+,hi,...,h a = 0, h^), where a € {1,...,N}. Then using the left 
action of the operator T> + (£ a — rj/2) we get: 

0(+) - n (_£ n/2) Slnh 11 TT C ° Sh 2C " 0) ~ C ° Sh 

+ [ ^ n ' J ^h(2e o + 7? )ll cosh 2 C W-cosh2C^ ) 

b^a 

x (+, hi, h a = 0, h N \+, hi, h a = 0, h N ) (4.15) 
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while using the decomposition (12.381) and the fact that: 



A + (-Ca + 7 ] /2)\+,h l ,...,h a = 0,...,h N ) =0 (4.16) 

it holds: 

V + (^-r ] /2)\+,h 1 ,...,h a = 0,...,h N ) = —^7 ^D + (-£a-r?/2) 

sinh(2£ a - t]) 

x \+,hi,...,h a = l,...,h H ), (4.17) 

and then we get: 

sinh 7] 



sinh(2£ a + t?) 

and so: 



D + {-^ a - V /2){+,h l: ...,h a = l,...,h N \+,hi,...,h a = l,...,h N ) (4.18) 



(+,h 1 ,...,h a = l,...,h N \+,h 1 ,...,h a = l,...,h N ) _ yr cosh 2(i 0) - cosh 2C^ h) 
(+,h 1 ,...,h a = 0,...,hu\+,hi,...,h a = 0,...,/i N ) Al _ CO sh 2cf 6) ' 



from which we have: 



(+,h!,...,h N \+,hi,...,h N ) _ -j-r 7/j 0) - vj 0) 

(+,0,...,0|+,0,..,0) H „(M_„N' 1 ' 

which proves the proposition, being 

(+,0,...,0|+,0,...,0) = n (0) 1 (0) . (4.21) 

l<6<a<N Va ~ Vb 

by our definition of the normalization N + . □ 



4.2 SOV-decomposition of the identity 

The following spectral decomposition of the identity I: 

2 N 

I = ^/i 4 |e,x- 1 (i))(e,^ 1 (f)|, (4.22) 
i=i 

can be given in terms of the left and right SOV-basis, where the fii = ((e, x _1 (i) \e, (i))) 1 is the analogous 

of the so-called Sklyanin's measur in our 6-vertex reflection algebra representations. Now using the result of the 
previous section we can explicitly write: 

l 

1 = E II ivi ha) -vi ha) )\^h 1 ,--- t h ti )(e,h 1 ,-'-,hu\- (4-23) 

ft,l,...,/i N =0 l<fe<a<N 

"Sklyanin's measure has been first introduced by Sklyanin in the quantum Toda chain 1691 . see also 1 120 1 and 1 121 1 for further discussions. 
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5 SOV characterization of T( A) -spectrum 



Here we construct in the £?±-SOV-representations the spectrum of the transfer matrix T(A) for the class of boundary 
conditions listed in Theorem 13. II Let us start giving the following characterization: 

Lemma 5.1. Let us denote with £7- the set of the eigenvalue functions of the transfer matrix T(X), then any r(A) € 
£7- is an even function of X of the form: 

t(A) = 2 sinh(A - rj/2) sinh(A + rj/2) coth C- coth (+ det M(m/2) 



+ 2 cosh(A - 77/2) cosh(A + rj/2) det M(0) 

9 

N 

+ sinh(2A - rj) sinh(2A + 77) 4 (cosh 2A) b ~ 1 . (5.1) 

6=1 

Proof. The transfer matrix 7~(A) is an even function of A so the same is true for the r(A) G £7- . Moreover, from 
the identities (13.291 ) and (13.521) after some simple computation the following identities are derived: 

T(±v/2) = a e (er//2) det M(0) = 2 cosh ri det M(0), (5.2) 
9 9 

T(±(r)/2 - ivr/2)) = iaJer)/2 - eivr/2) coth C+ det M(m/2) = -2 cosh 77 coth C- coth C+ det M(m/2), (5.3) 

q q 

both for e = +,—. This identities together with the known functional form of 7~(A) w.r.t. A imply the statement in 
the lemma. □ 

5.1 Transfer matrix spectrum in B -SOV-representations 

In this section we characterize the spectrum of the transfer matrix: 

T-(A) = r^ ) (A) + c + (A)/3_(A), (5.4) 

associated to the representations of the reflection algebra under the following class of boundary parameters: 

b + (A) = and 6_ (A) ^ 0. (5.5) 

Theorem 5.1. If the condition (|3,ip is satisfied, then 71(A) has simple spectrum and £71 coincides with the solu- 
tions of the discrete system of equations: 

T-(±Ci 0) )T-(±CP) = A_(C«)A4-d 0) ), Va e {1, N}, (5.6) 

in the class of functions of the form (|5.ip . where the coefficient A_(A) is defined by: 

A_(A) = a+(A)A_(A), (5.7) 

and satisfies the quantum determinant condition: 

a + (A)a + (-A + r?) detZY_(A - rj/2) = sinh(2A - r?)A_(A)A_(-A + rj). (5.8) 
9 
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I) The vector: 

1 N 

i r ->= e n^-^ M ) n (^ a) -^ 6) )i-^i,-^N), (5.9) 

h!,...,h N =Oa=l l<6<a<N 

defines, uniquely up to an overall normalization, the right T--eigenstate corresponding to t_(A) € E71. The 
coefficients in ()5.9[) are characterized by: 

Or_(d 1) )/Or_(d° ) ) = T_(d 0) )/A-(-d 0) ). (5-10) 

//) 77je covector 

1 N 

< r -i= e n^-«i M ) n (^ M -r/f 6) x-,/M,...^Ni, (5.H) 

/ii,...,/i N =0 a=l l<b<a<N 

defines, uniquely up to an overall normalization, the left T--eigenstate corresponding to t_(A) G £71. 77ie 
coefficients in ()5.11j) are characterized by: 

Qr.(CP)/QrM 0) ) = WM^J/Mtf'). (5-12) 

w/iere: 

qH = a + (d 1} ) = d + (d 1} ) = sinh(2g n + 2r ? ) i3) 
"d+(-d 0) ) a + (-d 0) ) sinh(2£ n - 2r/) ' 

Proof. In the £> -SOV representations the spectral problem for 71(A) is reduced to a discrete system of 2 N Baxter- 
like equations 

r_(cM* r _(h) =A_(e" ) )^_(T-(h)) + A_(-e"))^„(T r t(h)), (5.14) 
for any n G {1, N} and hG {0,l} N ,in the coefficients (wave-functions): 

V T _(h) = {-,h 1 ,...,h ti \T-), (5-15) 
of |t_) the 71-eigenstate associated to t_(A) G S71; here, we have used the notations: 

T^(h) = (/ii,...,/i n ±l,... J /i N ). (5.16) 

Being: 

A_(Cf)=A_(-C«) = 0, (5.17) 
the previous system of equations (I5.14p is equivalent to the following system of homogeneous equations: 



r_(d 0) ) -A_(-d 0) ) \ / * T -(h 1 ,...,h n = 0,...,h 1 ) \ = (0 

-a_(cP) r_(d 1} ) / 1 *r-(/ii, -,K = 1, -A) y lo 



(5.18) 



for any n G {1, N} with /i m ^ n G {0, 1}. The condition t_(A) G E71 implies that the determinants of the 2x2 
matrices in (|5.18p must be zero for any n G {1, N}, which is equivalent to (I5.6I ). Moreover, the rank of the 
matrices in (I5.18P is 1 being 

A-(-d 0) )/0 and A_(C«)^0, (5.19) 
and then (up to an overall normalization) the solution is unique: 

V T _(hl,...,h n = l,...,h 1 ) _ T_(Cf) 



V T -(h u ...,h n = 0,...,h 1 ) A_(-d' 



(0 V 
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for any n G {1, N} with h mj t n G {0, 1}. So fixed t_(A) £ there exists (up to normalization) one and 
only one corresponding 71-eigenstate |r_) with coefficients of the factorized form given in (|5.9[) - ()5. lQj) : i.e. the 
71 -spectrum is simple. 

Vice versa, if r_(A) is in the set of functions (15.11) and satisfies (15.61) . then the state |r_) defined by (|5.9p - (|5.10p 
satisfies: 

(-, h x , h N | T-(Ci hn) )\r-) = r_(d M )<- h, h N \r-) Vn G {1, N} (5.21) 
for any £>_-eigenstate (— , hi, /in I and this implies: 

(-,h 1 ,~M\T-(\)\T-)=T-(\){-,h 1 ,~'MT-) VAGC, (5.22) 
i.e. t_(A) G S71 and |r_) is the corresponding 71-eigenstate. 

Concerning the left 71-eigenstates the proof is done as above we have just to remark that in this case the matrix 
elements: 

<T-|T-(d M )|-A,->N), (5.23) 
can be computed by using the right B- -representation: 

r_(Ci M )^_(h) =D_(d M )^_(T;(h))+D_(-e" ) )^ r _(T+(h)), VnG{l,...,N} (5.24) 

where: 

V t _Ql) = (t-\-M,~.M, (5-25) 

and the coefficient D_(A) reads: 

D_(A) = d+(A)D_(A). (5.26) 

□ 

It is worth pointing out that for the analysis of the continuum limit it is interesting to get a reformulation of this 
characterization by functional equations. The construction of a Baxter Q-operator can play an important role to 
achieve this aim. Let us recall that a Q-operator is of a one-parameter operator family which satisfies properties of 
the type: 

[71(A), Q(A)] = 0, [Q(X), Q( M )] = 0, 71(A)Q(A) = a (A) Q{X/q) + /3 (A) Q(Xq) , (5.27) 

where a (A) and /3 (A) are some characteristic functions of the constructed Q-operator. Indeed, if this functional 
equation coincides with the discrete system (|5.14p in the spectrum of the £>±-zeros we can use the Q-operator to 
reformulate by functional equations the SOV spectrum characterization. A Baxter Q-operator has been constructed 
irl^l 1 621 l63l and ll64l for some classes of representations of the reflection algebra; then, it will be important to make 
a connection with the present SOV analysis. Let us observe that for root of unit rj = 2mp'/p (p and p' G Z-°) the 
existence of non-trivial solutions of the Baxter equation leads to the functional equation: 

det£>(A) = 0, AgC (5.28) 

v 

20 See also 1941 and 1 124 1 for the construction of the Q-operator in the spin 1/2 and higher spin XXZ quantum chain with twisted boundary 
conditions. 
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where D(\) is the following p x p matrix: 

/ r_(A) -/3(A) 
-a(q\) T-(q\) 



D(X) = 





-f3(qX) 



{-Piq^X) 



a(qP~ 2 X) T^(qP- 2 X) 











-a(A) \ 




aiqP^X) r_(gf- x A) / 





-P{qP- 2 X) 
p— 1 



(5.29) 



written only in terms of the T- -eigenvalue t_(A) and its determinant is clearly a function of A = A p . Note that the 
method based on the combined use of the fusion of transfer matrices 111251 11261 and the truncation identity for root 
of unit rj, lfT27l and |[T28l[T29l . should lead to the same equation dOH]) . 



5.2 Transfer matrix spectrum in # + -SOV-representations 

In this section we characterize the spectrum of the transfer matrix: 

T + (A)^ (+) (A) + c_(A)S + (A), 



(5.30) 



associated to the representations of the reflection algebra under the following class of boundary parameters: 

b_ (A) = and b + (A) ^ 0. (5.31) 
Let us now write the left and right eigenstates of the transfer matrices (I5.301 l in the £> + -SOV-representation: 

Theorem 5.2. If (|3.ip is satisfied, then 7+ (A) has simple spectrum and coincides with the solutions of the 
discrete system of equations: 

r + (±d 0) )r + (±Ci 1) ) = D + (-tf>)D + (<f), Va e {1, N}, (5.32) 

in the class of functions of the form (15. 1 D , where the coefficient D + (A) is defined by: 

D+(A) = d_(A)D+(A), (5.33) 

and satisfies the quantum determinant condition: 

d_(A - rj/2)d_(-A - rj/2) detU + (X) = sinh(2A + 2t/)d + (A - r//2)D+(-A - r?/2). (5.34) 



(5.35) 



I) The vector: 

1 N 

|r+)= £ UQrM ha) ) II (vi ha) -V ( b hb) )\+,h 1 ,...,h N ), 

hi,...,h N =0a=l l<6<a<N 

defines, uniquely up to an overall normalization, the right T+-eigenstate corresponding to t+ (A) € . 77ie 
coefficients in (|5.35f) are characterized by: 



Or + (Ca W )/Or + (tf ) )=r + (tf))/D + (CW)). 



(5.36) 
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II) The covector: 

1 N 

(r+\= £ U&M ka) ) II (^-r/^X+^i,...,^!, (5.37) 

hi,...,h N =Oa=l l<6<a<N 

defines, uniquely up to an overall normalization, the left T+-eigenstate corresponding to r + (A) € £7-+- The 
coefficients in ([5.37P are characterized by: 

Qr + (Ci iy )/QrM 0) ) = r + (Ci 0) )/ (aW4 +) D + (-Ci 1) )) , (5-38) 



(+) _ a_(d 0) ) , d_(d 0) ) sinh(2£ n - 2g) 



„(+) - Jz(<f 

d_(-cP) a_(-CP) A;i +) sinh(2£ n + 2r/) 

Proof. Taken a 7+-eigenstate |r + ) corresponding to the eigenvalue r+(A) € E71, it has coefficients: 



V T+ (h) = (+,h 1 ,...,h N \T+), (5.40) 
in the SOV-basis, which satisfies the following discrete system of Baxter-like equations: 

r + (C^ T+ (h) =D + (C^))* T+ (T+(h)) + D + (-C^)* r+ (T-(h)), (5.41) 
for any n £ {1, N} and he {0, 1} N . We can rewrite this as the following system of homogeneous equations: 



r+(d 0) ) -D+(d 0) ) \ / * T+ (h 1 ,...,h n = 0,...,h 1 ) \(0 



(5-42) 



being: 



D+(-d 0) ) = D+(d 1} ) = 0, D+(C<°>) ^ and D+(-tf>) 0. (5.43) 
From which the condition t+(A) € 2-7-+ directly implies (15.61 ) and moreover it has to hold: 

* T+ (h u ...,h n = i,...,h 1 ) _ r + (d 0) ) 



y T+ {h 1 ,...,h n = o,...,h 1 ) D+ (ci 0) )' 



(5.44) 



for any re € {1, N} with h m + n G {0, 1}. This fixes the factorized form given in (|5.35|) - (|5.36p for the T+- 
eigenstate |r + ) and implies the simplicity of the 7+-spectrum. Taken a r+(A) solution of (15.61 ) in the class of 
function (15.11) and constructed |r+) by (|5.35p - ()5.36p . then the proof that t + (A) € £7^. and |r + ) is the corresponding 
7+-eigenstate can be given following the same steps presented in Theorem 15. II Concerning the left 7+-eigenstates 
the construction is done as above we have just to remark that in this case the matrix elements: 

(r + \T + (Cj l h - ) )\+,h 1 ,...,h N ), (5.45) 

can be computed by using the right B + -representation: 

r+(C^)f T+ (h) =A + (C^ ) )^ + (T+(h))+A + (-C^ ) )^ T+ (T-(h)), Vne{l,...,N} (5.46) 

where: 

9 T+ (h) = (T+\-,h 1 ,...,hn), (5-47) 

and the coefficient A + (A) reads: 

A+(A) = a_(A)A+(A). (5.48) 

□ 
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6 Scalar Products 



The presentation will be done simultaneously for £> e -SOV-representations with e = +, — . 

Proposition 6.1. Let (a e \ an arbitrary covector and let \/3 e ) be an arbitrary vector of separate forms: 

1 N 

(«e|= E II^^ M ) II (vi ha) -vi hb) )(e,hi,..,hu\, (6-1) 

hi,...,/i N =0 a=l l<6<a<N 
1 N 

IA>= E II^^) II (^ M -r 7 ? i ' ) )|6,/ il ,...,/ tN ), (6.2) 

hi,...,/i N =0a=l l<6<a<N 

/n ?/ze B t -eigenbasis, then the action of (a e \ on \/3 e ) reads: 

(a t m = det ||A<J A) || wiA s^^-C^M^)^)^- ( 6 - 3 ) 



h=0 



Moreover, ifr e (X) ^ T~' e (\) € X-r E , f/je action of the %-eigencovector (r e | o« ?/ze 7^- eigenvector \r' t ) is zero, and in 
particular it holds: 

N 

E^S'^' 7 "^ Va€{l,...,N}, (6.4) 



(r t') 

where the c b £ ' are defined by: 



N 

r e (A) - r £ '(A) = sinh(2A - r?) sinh(2A + 77) E 4 T£ ' T£,) ( cosh 2A) 6 " 1 . (6.5) 



6=1 

Proof. The formula (j4.7|) and the SOV-decomposition of the states (a € \ and \/3 e ) imples: 

N 

II ) ■••) '/|\| ; 

/ii,...,/iN=0 a=l 



(aei^) = e n^ ,-,^ 5 )!!^^)^^), (6-6) 



where ...,xn) = ni<6<a<N( x a ~~ x &) * s tne Vandermonde determinant which for the multilinearity of the 

determinant implies (|6,3p . 

The 7^-eigenstates (r e | and [r e '} are left and right separate states then the action (t 6 \t^) is also given by (|6,3p and to 
prove (t £ |t € ') = for r e (A) 7^ t € '(A) € we have just to prove (I6.4p . It is simple to remark that: 

°' & ^ fro smh(2^ ) - ?? )sinh(2^ + r ?) ' ^ 

so in the case e = — , we can use the equations (15.101 ) and (15.121 ) to rewrite: 

-A_(c«)Q ri (ci 0) )g._(d 1) ), (6-8) 

and 

Qrl(d 0) )Qr_(d 0) )(T-(d 0) ) " Ti(d 0) )) = (^aa)" 1 [A-(CP)Qri (C^Qr. (CP) 

- fc a a a A„(-C(°))Q Ti (C( 1) )0r_(d 0) )]. (6-9) 
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Then by substituting them in (16.71) we get (I6.4I ). Similarly, in the case e = +, we can use the equations (15.361 ) and 
(15.381 ) to rewrite: 



- 4+)ai + )D + (-C( 1 ))Q T ;(d 0) )Q. + (d 1) )], (6-10) 

and 

Qt' + (Ci 0) ) Qr+ (Ci 0) ) (t+ (Ci 0) ) - r' + (C o 0) ) ) = ^D^jg,,^)^^)) 

-D + (c(°))Q < (d 1) )Q. + (d 0) ). (6.1D 



Then by substituting them in (16.71 ) we get (I6.4I ). □ 

It is worth remarking that the vector ((e, a\y G is of separate form in the Q-eigenbasis thanks to the Hermitian 
conjugation properties of the Yang-Baxter generators. Then the previous result describes also the scalar products for 
these states. The determinant formulae obtained here can then be considered as the SOV analogous of the Slavnov's 
scalar product formula H122IU23l - ll34l which holds in the framework of the algebraic Bethe ansatz. 



7 Boundary reconstructions of strings of local operators 

Here, we present identities between couples of (boundary) generators of the reflection algebra and (bulk) generators 
of the Yang-Baxter algebra which can be used to reconstruct any local operator in terms of the boundary operators. 
This program is explained and developed for simple strings of local operators likes: 

7.1 Mixed bulk and U± -boundary reconstructions 

Let us start recalling the bulk reconstruction formulae: 

Proposition 7.1 ([34]). Let x n £End(R n ) be the generic local operator in the local quantum space R n , then it 
admits the following reconstruction in terms of the generators of the Yang-Baxter algebra: 

n-l Tfr^h n_1 

x n - II J tfi ^ o(M (Cf )gfc?og) II r-W) (7-2) 

N T(7^ N 

a=n+l ' vS / a=n _|_i 

= II T (^ 1) )^o(M (C«)xo)-^| 7 f- I] r- 1 ^) (7-4) 

a=l g lUJ a=l 

= n ^^^ det M(n tro(M ° (c " ))a;o) n t <&\ ™ 

a=n+l " \^> n / a=n _\_i 

where we have used the identities: 

T((^)T(C^) = detM(U). (7.6) 
q 
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In II 1301 a reconstruction of local operators which use the reconstruction of the propagator by the (bulk) transfer 
matrix T((i ha ') of the Yang-Baxter algebra and the elements of the (boundary) reflection algebra U + (A) has been 
derived. In this subsection we reproduce this result and provide other three equivalent reconstructions: 

Proposition 7.2. Let x n £ End{R n ) be the generic local operator in the local quantum space R n , then it admits the 
following reconstructions in terms of the generators of the (boundary) reflection algebra U+ (A): 

n-l =. ( JQ)s n-1 

n-l =- (Al)\ n-l 

a=i det q U + (U) ^ 

and the following ones in terms of the generators of the (boundary) reflection algebra U— (A): 

N =- , (1). N 

a=n+l 9 a=n+l 

N =. ,J0), N 

a=n+l 9 \^ n / a=n+l 

//ere, we /jave denoted^: 



7±(A) = MA)At(A) + M-AM±(-A), (7.1 1) 



ant/ 



a±(A) = a ± (A)| c+=W2 = d±(A)| c+=i7r/2 = cosh(A t ty/2) ^ , (7.12) 

detZ7_(A) = A_| c+=i7r/2 (A + 7 ? /2)A_| c+=i7r/2 (-A + 7 ? /2), (7.13) 
detW + (A) = D + | c _ =W2 (-A-? ? /2)D + | f _ =W2 (A-r ? /2). (7.14) 

Proof. The first reconstruction in terms of U + (A) is the result proven in Proposition 1 of 111301 : similarly it is 
possible to prove our second reconstruction in terms of U + (A). Let us prove here the reconstructions in terms of 
U- (A); by definition of propagator operator it holds: 



[] T(cP)tr {U^(X)cjlx^al) = (-l) N tr (L ,„(A) . . . L i(A)L 0N (A) . . . L 0n+1 (\)K_(\)a y 

0=1 

n 

2$ n+1 (-A) . . . 2# N (-A)L& (-A) . . . 4U-A)4'^) [J nC«) (7-15) 



0=1 



where Loa(A) — ^Oa(A — Ca — f° r an Y a € {1, N}. Now using the identities: 

L n(d 1} ) = Pon sinhr?, L n(d 0) ) = - smh V a y Pfaal (7.16) 



T-(A) = tr (U-(X)K + (\)\ (+=iv/2 ) = cosh(\ + r)/2)(A-(\)+V-(\)), 
T+(A) = ^ (W + (A)if_(A)| c _ =W2 ) = cosh(A-r,/2)(^ + (A) + O + (A)). 
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where Po n is the permutation operator between the 2-dimensional spaces Ro and R n . Then, the r.h.s of (T7.15t 

n 

tr (U-(-(^))x n llT((V), (7.17) 



computed in A = — Cn reads 



a=l 

being: 

4UCP>r>0 = sinhT? [x P 0n ] t0 a y = S mh V P^a y x n . (7.18) 
Similarly, the r.h.s of (17.151 ) computed in A = — (4 reads: 



XntroiU-i-C^YlTi&y), (7.19) 

a=l 

being: 

4°^o^On(d 0) ) = - sinh?? [Ponxof o% = x n a v L 0n ((W) . (7.20) 

By using these formulae the reconstructions (|7.9t and (17.101 ) simply follows. □ 

The previous proposition naturally implies the following: 

Corollary 7.1. The following annihilations identities hold: 
I) For the generators of the reflection algebra U- (A): 

A_(Cf )C_(±d 1) ) = A-(Ck 0) )^-(-d 1) ) = 0, (7.21) 

A_(-CP)C_(±C^) = A-C-tfJ)©-^ ) = 0, (7.22) 

2?_(Cf )B-(±C«) = P_(C(° ) M-(-CP) = 0, (7.23) 

P_(-C«)5_(±Cf) = P_(-d 1) M-(d° ) ) = 0, (7.24) 



ant/ 



s_(±4 0) )^-(±c^) = m+c^V-k^) = £-(±d 1} M-(d 0) ) = o, (7.25) 

C_(±d 0) )C-(±d 1} ) = C_(±d 0) = C_(±C«)P-(C( 0) ) = 0. (7.26) 



//) For £/2e generators of the reflection algebra U+ (A): 



^ + (-Cf)S + (±Ci 0) ) = ^ + (-a O) )^+(C«) = 0, (7.27) 

AfCCW^+Ci^) = ^l + (d 1) )P + (-d 0) ) = 0, (7.28) 

V + {-(g>)C+{±<W) = P + (-CfM + (d 1) ) = 0, (7.29) 

P + (C«)C + (±Cf) = P + (C( 1) )^ + (-Cf) = 0, (7.30) 

B+(±Cf )B + (±C«) = B+(±<f )X> + (d 1} ) = S + (±d 1) )^+(-Cf) = 0, (7.31) 

C + (±Cf )C + (±CP) = C + (±C£V+(C«) = C^itf ))Ah(-<< >) = 0. (7.32) 
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Proof. The proof of these annihilation identities can be done along the same line used for the bulk case. Let us 
sketch it; by using the reconstruction formulae of the previous proposition and the identities: 

T± (Ci 1} )T± (Cf ) = det U± (f „ ) , (7.33) 



it also holds: 



nrwnn<P) - ^M&^p^ (7 , 4) 

a=i a=i det,W+(eO 

r+riVo^+t-cf )4 1) )^o(^+(d 1) )4 2) )r+(d 0) ) 



(det W+(£ n )) 2 



(7.35) 



and: 



TT ntf 1 )^ TT r-^cP) = tr °^-tf ) )^ ) )^f-(-^ 1) )^ ) ) (7.36) 

11 VSa ii VSa ; det W__(£ n ) 



a=n+l a=n+l 

_ T- )*r (M-C^rf% (M_ (Cf rf } )T- (Cf ) 



(det 9 W_(^)) 2 



(7.37) 



where we have defined: 



4 3) - y W = at (y«y«f °l (7-38) 



we can use now these formulae to derive all the annihilation formulae. □ 
7.2 U± -boundary reconstruction 

Let us remark that the bulk reconstruction formulae of Proposition 17.11 the corresponding annihilation identitieJ^l 
and the Yang-Baxter commutation relation implies that the generic string of local operators of the form: 

n 

J] (7-39) 

a=l 

can indeed be represented as linear combinations of the following strings of bulk operators: 

(7.40) 

then the following identities between bulk and U + -boundary generators: 



r r(A°)\ q- try 

tr (M ((W)y ) - -^-i . = tr (U + ($>)yQ) j—hi r , (7-41) 



-(OK -7- c^(O) 

detgM(4„) det £Y+(£ n , 
imply that we can also write (17.391 ) by a linear combinations of the following strings of U + -boundary generators: 



(7.42) 



22 See Lemma 5.1 of 1291 . 
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Similarly, the identities between bulk and Z//_ -boundary generators: 



MMo(d%D) ,f^ = ir o(W-(C< 0, )l/o)^T (7.43) 
det 3 M(£ n ) det q U-{£ n ) 



imply that the generic string of local operators of the form: 

N 



H x a (7.44) 



a=n 



can indeed be represented as linear combinations of the following strings of U- -boundary generators: 

(7-45) 

Here we show explicitly as these reconstructions work for some special string of local operators: 

Proposition 7.3. Let us consider the open XXZ spin chain with general K- and diagonal or triangular K + , then 
the following boundary reconstruction holds: 

( ^N+i-n TT MCP) tt sinh(£ a + & - rj) 
n ' ' ' n = ("I) llr^OK 11 sinhfe + &) 

a=n J n<a<b<N VS S0/ 

x^ ( ci 0) ) • • • g-(ff ? ) :rr^ • • • JTrnk - « M > 

det ? W_ (£n ) detq U- (£ n ) 

Lef consider the open XXZ spin chain with general K + and diagonal or triangular K^, then the following 
boundary reconstruction holds: 

A <Md 0) ) t-t stahfo + & + 



1 " fid-tf^^ 

(0), r+(c (o ); 



: b + (cJ 1) )--^ + (C«)-t5# 



det g ZY+(£ n ) det 5 W+(6)' 



(7.47) 



Proof. Let us prove explicitly the first reconstruction as for the second one can proceed similarly. First of all by 
using the bulk reconstruction of Proposition 17. 1 l one gets: 

^- *n = (-l) N+1 -^(Ci 0) ) • • • B(C^) f ( 5f ) • • • (7.48) 

detqMo(£ N ) det g M (4„) 

Now by using the bulk annihilation identities B((^ )D(^') = 0, we get the set of identities: 

B(CS } ) ■ ■ ■ B(C<!»)T«£>) = B«< ») ■ ■ ■ B(^>)A(C^) 

= n ^^Ftt 4 ^'^ 3 ^' • • • B(cf,) 

= ff ^^z^r ^^mcg^B^) • • • ^(d o) ) (7.49) 
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where to commute and ) we have used the Yang-Baxter commutation relation: 

BK f )^(c< 1 ') = **&ftBffL M <Pm<m) + -^L.ldVKi") (7.50) 

sinh(£ n - £ N ) smh(£ n - £ N ) 

which with the bulk annihilation identity 5(^4 )-S(Cn ) = imply the above second identity, the third one is 
obtained by reiterating the commutation and finally the forth one by using once again the annihilation identity 
B((^)D((^) = 0. Repeating the same procedure to commute the others transfer matrices through the product of 
the /^-operators we get: 

a „- = f_nN + l-n TT smh^-j^) Bj^Tj^) g((gWg) 

N n 1 J sinli(£ a -£ 6 ) det g M (£ N ) det ? M (£„)' 

and then by using the boundary -bulk identities (|7.43p we can rewrite it in the form: 

„- n -- ( _^ + i-n tt sinh(g a - 6 - ry) g-ccf^cg) g-(d 0) )r-(d 1} ) r7 , 9 . 



N 



, N sinh(£ a -£ 6 ) det ? W_(£ N ) det 9 W_(£ n ) 



We can now use the boundary annihilation identities (|7.25p and the reflection algebra commutation relations to move 
all the transfer matrices 71 (Ca ) to the right and transform them into 71 (Ca )• More in details, from the annihilation 
identities and the definition of the transfer matrices it holds: 

BAC^T-m = ^^B-^T-m, (7-53) 

a +(Cn ) 

and so the first transfer matrix on the right in (|7.52p can be rewritten in the desired form. Now let us apply the 
procedure to the product B_(^j7L(c2i) B -(Cn 0) ) to rewrite it in the desired form B_(C^i)S-(Cn 0) )7l(ci+i), 
the annihilation identity (|7.25p implies: 

B-«Si)T-«2i)B-(^ 0) ) = «+(CSi)^-(CSi)^-(CSi)e-(a o) ) (7.54) 

_ />(l) v sinh(g n - £ n+i ) sinh(g n + £ n+1 - rj) 
° +lWlJ smh(£ n - Cn+i - V) sinh(^ n + U+i) 
xB^l)B-(Ck 0) )A-(&i)- (7-55) 
In the second identity we have used the reflection algebra commutation relation: 

.4-(e)s-(e>) - = tf" z t s "* l) i S ""1; + !T 7' «-«. m M-«&) 

smh(£ n - Cn+i ~ V) smh(£ n + U+i) 



sinh(2£ n — 2rj) sinh 77 



sinh(£ n - £ n+ i - 77) sinh(2£ n - rj) 



- ■v^TVhax ^W) (7 " 56) 

smh((„ + £ n+ i) smh 2Ai 

and the fact that due to the presence of the second and third terms on the right of this formula give zero 

when inserted in (|7.52|) . Now by using the commutativity of the £>_ -generators and once again the annihilation 
identities it holds: 

Md2i)*-& 0) )^(<&) = B -^ B -^ T ~^\ (7.57) 

a+(Q ) 

so that we have accomplished our task for N = n + 1 having the product 7l(Q+i)7l(Cn ) to the right of (|7.52p . 
Reiterating this procedure for the remaining transfer matrices 71 (Ca^) for a € {n + 2, N} we obtain our result. 

□ 
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8 Matrix elements 

8.1 From B -SOV representation 



Here we consider the transfer matrices (15.41 . then the following proposition holds: 

Proposition 8.1. Let (r_ | and \t'_) be a generic couple of left and right T--eigenstates, then we have: 

(K—e T ~ sinh r\ \ ^ n 
_ 2" sinh C- ) -TT Sinh(^ a + &-??) ,, (-n,T_,ri) 

Vi-|«„---« N r'-/- (o) (0) 11 s inh(C a + 6) 2N-n l|2ja ' 6 11 

v \J]n J n<a<b<N VSa S0; 

A. 0+(d 1) )7i(d 1) )Qr_(d 1) )Qr' (Ci 1) )smh2^ a 

x n m — — , (8.i) 

where ||S Q fe ' n ' T ' T ~^|j is the (2N — n) x (2N — n) matrix of elements: 

forae {l,...,n-l}, be {1,...,2N - n}, (8.2) 
} for a G {re, N}, b G {1, 2N — n}, (8.3) 
} /or a G {N + 1, 2N - n}, be {1, 2N —n}. (8.4) 

Proof. Here we use the reconstruction (I7.461 l for cr~ • • • cr^ , then we can act with the product of transfer matrices on 
the right 7^-eigenstate \t'_) and we are left with the following computation: 

S_(4 0) )---£_(Cf )|Ti>. (8.5) 
From the decomposition of \t'_) in the ,B_-eigenstates, we get: 

N 1 n-l N 

s-(ci o) ) • • • b-«< o) )k-> = n >) e n ^(cf a) ) n b „„ ,, ^ 

a=n hi,...,ft„_i=0 o=l o=n 

V( V [ h l} , .., ^r l} , 77W, r?^)!^, ^-i, 1, 1). (8.6) 
Let us rewrite the £>_ -eigenvalues in terms of the 77^ : 

b-,(a> s * 2 f s :r ( h(2A - n H« - ■ » 7 > 

a=l 



(-,n,T_ 










y.) 




r 


V (-.™: T - 


S-) 







and then we have: 



n B -^ 1 ,..^,..,i}(ci o) )^rf i v..,eT i) ) ^ i) ,--- ) tf ) ) 



/ 1 \N — n T ■ 1. \ N— 71 

( — 1) K—e — sinh7; \ ^ 



2 N sinhf- 



]J sinh 2(6,-7?) 



x^V.erUf,.,^,^,.,^). (8-8) 
Using this last formula and taking the scalar product we obtain our result. □ 
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8.2 From /3+-SOV representation 

Here we consider the transfer matrices (|5.30p . then the following proposition holds: 

Proposition 8.2. Let us consider the open XXZ with transfer matrix 7+ (A) (|5.30|) and let (r+| and\T' + ) be a generic 
couple of left and right T+-eigenstates, then we have: 



I (-l) N K + e r + sinhr? 

" j_(d 0) )T+(d 0) )<9r + (d 0) )Qri(d 0) )sinh2e a 

x n ^ ^ , (8.9) 

where ||S + 1 1 « f/je (N+n) x (N+n) matrix of elements: 

E^' ft,T+y+) = (7 7 «) (6 " 1) / rae{l,..,n}, 6 G {1, N + n}, (8.10) 

E (+,n,r+y + ) _ / (o^y 6 - 1 ) /oraG{n + i,...,2n}, 6 G {!,..., N+n}, (8.11) 



s (+.»,r+,r;) ^ M <r+£+) for a e {2n + 1, N + n}, 6 G {1, N + n}. (8.12) 

Proof. Here we use the reconstruction (T7.47t for cr^ • * • a~, then we can act with the product of transfer matrices on 
the right 7+-eigenstate [r^_) and we are left with the following computation: 

B + (C? ) )---B+(^ 0) )I<>. (8-13) 
From the decomposition of \t' + ) in the ,B + -eigenstates, we get: 

n IN n 

b+({W) ■ ■ ■ b + ((U)\t> + ) = n^;(cP) e n ^ + (^ ) )iiB +l{ o 1 ... 1 o,fc^ 1 ,..,fc N }(d i) ) 

a=l /i n _|_i,...,/ifj=0 a=n+l a=l 

V(77< 0) , ...^^^Shi 1 - , -..^^lO, ...,0, fe N >. (8.14) 

Let us rewrite the B + -eigenvalues in terms of the 77^ : 



2 N sinhC- 

a=l 



and then we have: 



(-l) N K + e T + sinh?? 



2 N sinh £ 4 



— ^- Yl sinh 2£ a sinh2(£ a + n) 



xn1^ 1 ^■.•^ l ^^ ^.■^^nto l) .■■■.'A* N) )■ ( § - 16 ) 

Using this last formula and taking the scalar product we obtain our result. □ 
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9 Conclusion and outlook 



We have analyzed the integrable quantum models associated to the transfer matrices corresponding to one general 
non-diagonal and one diagonal or triangular boundary matrices. For these integrable quantum models, defining the 
open spin 1/2 XXZ quantum chain in the same class of non-diagonal boundary matrices for the homogeneous limit, 
we have obtained the complete SOV-characterization of the transfer matrix eigenvalues and eigenstates, the proof of 
the simplicity of the spectrum and determinant formulae of N x N matrices for the scalar products of separate states. 
Finally, matrix elements of a class of quasi-local operators have been computed on the transfer matrix eigenstates in 
determinant form by using the reconstruction of these operators by the Sklyanin's quantum separate variables. The 
relevance of these findings in the framework of the non-equilibrium systems like the partial asymmetric simple ex- 
clusion processes (PASEP) will be described in 111311 and for the most general symmetric simple exclusion processes 
in H132II . where moreover further matrix elements of quasi-local operators will be computed. 

In the literature of quantum integrable models there exist different applications of separation of variable methods for 
computing the matrix elements of local operators. An important example is presented in the Smirnov's paper [ 120], 
where determinant formulae for the matrix elements of a conjectured basis of local operators have been derived in 
Sklyanin's SOV framework for the quantum integrable Toda chain [69]. There is a strong analogy among Smirnov's 
formulae, those that we have here derived and more in general those which appear in the series of papers [|2|, ||84] l85l 
and IIIHS. The main differences in all these formulae are due to model dependent features, like the nature of the 
spectrum of the quantum separate variables. In fact, it is worth citing also the results of the papers M_33. 134] on 
the form factors of the restricted sine-Gordon at the reflectionless points in the S-matrix formulatiorr 3 l The form 
factors there derived^! can be represented once again as determinants and the connection with SOV emerges on the 
basis of the semi-classical analysis of 11 3311 - there also used as a tool to overcome the problenj^l of the local fields 
identification. 

Let us comment that in this paper we have followed a different approach for the reconstruction of local operators 
w.r.t that used in ||29]|3TI . A part the different framework, ABA in |[29l [3TTl and SOV in this paper, we have decided 
to reconstruct local operators directly by the quantum separate variables of the 6-vertex reflection algebra and not 
in terms of those of the 6-vertex Yang-Baxter algebra. The main motivation to do so is related to the increased 
complexity of the functional relations among the generators of these two algebras in the general non-diagonal cases 
which make more complicated compute the action of the quantum separate variables of the 6-vertex Yang-Baxter 
algebra on the eigenstates of the 6-vertex reflection algebra transfer matrices. Our current approach is of course 
more natural as the action of the quantum separate variables of the reflection algebra on the corresponding transfer 
matrix eigenstates has a simpler form. However, it is worth commenting that this reconstruction program is not yet 
completed as we have so far constructed explicitly only some classes of quasi-local operators by our approach but we 
believe possible to use this type of reconstruction to compute all matrix elements of local operators. One important 
motivation is to derive form factors of local operators expressed by determinant formulae as it was obtained in ||2| 
for the 6-vertex transfer matrix with antiperiodic boundary conditions. Indeed, the knowledge of the form factors 
of local operators is an important step toward the complete solution of the quantum model as the form factors 
represent an efficient numerical tools for the computation of two point correlation functions. In fact, we can rewrite 



23 See 1 135|-| 139| and references therein. 

24 Note that recently in [ 140] these results have been connected to the important achievements obtained in 1 140 ]-[ 148] where a fermionic 
basis of quasi-local operators has been introduced in the infinite volume limit of the XXZ spin 1/2 chain. 

25 Let us recall that this is a longstanding problem in the S-matrix formulation. The description of massive IQFTs as (superrenormalizable) 
perturbations of conformal field theories 01491 - 111531 by relevant local fields 1 1541- 11571 has been at the origin of the attempt of classifying the 
local field content of massive theories (the set of the solutions to the form factor equations [ 158. 159]) by that of the corresponding ultraviolet 
conformal field theories. Several results are known which confirm this characterization, see for example [ 160 1-[ 163 | and the series of works 
fl64l-IT67l. 
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correlation functions in spectral series of form factors and then we can try to use the same approach developped in 
111681 in the ABA framework and used in the series of worksl^l II176II - II181II . This is a concrete project as also in our 
SOV framework will be possible to have representations for the scalar products and complete characterization of the 
transfer matrix spectrum in terms of solutions of a system of Bethe equations. 

Finally, let us comment that the analysis developed in this paper define the required setup to extend the results on the 
spectrum characterization and the scalar product formulae in the SOV framework to the most general non-diagonal 
spin- 1/2 open XXZ and XYZ quantum chains. Indeed, the so-called gauge transformations! 27 ! can be used also in 
the reflection algebra framework to reduce the spectral problem to one analyzable by SOV. More in details, both the 
transfer matrices of 8-vertex and 6-vertex reflection algebras associated to the most general integrable boundaries 
matrices can be reduced by gauge transformations to those of a dynamical 6-vertex reflection algebra of elliptic and 
trigonometric type, respectively, with one triangular boundary matrix. The implementation of the SOV analysis for 
the spectral problem of these dynamical 6-vertex systems is currently under study in collaboration with N. Kitanine 
and it consists in the generalization to the dynamical case of the SOV results derived in this paper for the standard 
reflection algebra. It is then worth mentioning that in [85 1 85 ] it will be shown as the SOV results from the spectrum 
up to the form factors of local operators can be extended from the standard Yang-Baxter algebra to the dynamical 
one. 
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